Approaching ideal weak link behavior with three dimensional aluminum nanobridges 
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We present transport measurements of unshunted do superconducting quantum interference de- 
vices (SQUIDs) consisting of 30 nm wide aluminum nanobridges of varying length L contacted with 
two and three dimensional banks. 3D nanobridge SQUIDs with L < 150 nm (approximately 3-4 
times the superconducting coherence length) exhibit ~ 70% critical current modulation with applied 
magnetic field, approaching the theoretical limit for an ideal short metallic weak link. In contrast, 
2D nanobridge SQUIDs exhibit significantly lower critical current modulation. This enhanced non- 
linearity makes 3D nanobridge Josephson junctions well suited to optimize sensitivity in weak link 
SQUID magnetometers as well as realize ultra low-noise amplifiers and qubits. 



Josephson junction circuits based on weak links are 
well suited for applications requiring high transparency, 
small area junctions. For example, superconduct- 
ing quantum interference devices (SQUIDs) constructed 
from nanoscale bridges are widely used in magnetiza- 
tion measurements of magnetic crystals and ferromag- 
netic clustersi^— . Furthermore, weak link junctions are 
attractive candidates for coherent quantum circuits since, 
by construction, they are free of lossy tunnel oxide — a 
potential source of decoherenceA. However, these perfor- 
mance gains come at the expense of reduced nonlinearity 
in the current-phase relation (CPR) when compared to 
conventional tunnel junctions. The CPR in weak link 
devices depends strongly on the junction geometry^. In 
particular, thin film bridges with lateral dimensions com- 
parable to the superconducting coherence length ^ are 
predicted to have maximal nonlinearity when contacted 
with banks of greater thickness than the bridge (3D), as 
opposed to banks of equal thickness (2D). The 3D banks 
enhance nonlinearity by acting as good phase reservoirs 
and confining the variation of the superconducting order 
parameter to the weak link§. 

In this letter, we present measurements of the mag- 
netic flux response of unshunted dc SQUIDs comprised 
of aluminum nanobridges connected to 2D and 3D banks. 
Our measurements show enhanced critical current mod- 
ulation in 3D nanoSQUIDs and confirm that thin film 
aluminum nanobridges up to 150 nm in length (w 3 — 4^) 
with thick contacts approach the nonlinear CPR of an 
ideal, short, metallic weak link. Longer devices and those 
with 2D banks exhibit reduced modulation, in quantita- 
tive agreement with our numerically computed CPRs. 
Thus, 3D nanobridge junctions are robust nonlinear ele- 
ments for use in nanoSQUID magnetometers, amplifiers, 
and coherent circuits such as quantum bits. 

Our devices consist of a small superconducting loop 
interrupted by two nanobridges. Each bridge is an 8 nm 
thick, 30 nm wide aluminum film of variable length from 
75-400 nm. These dimensions can be reliably achieved us- 
ing electron beam lithography and a liftoff process. For 
2D devices, a single electron-beam evaporation is per- 
formed at normal incidence through a bilayer shadow 
mask to simultaneously define a narrow bridge and wide 
banks of the same thickness. To create 3D banks, we tilt 




FIG. 1: (Color online) Atomic force micrographs of 8 nm 
thick, 30 nm wide aluminum nanobridges contacted with (a) 
8 nm thick 2D banks and (b) 80 nm thick 3D banks. Two 
nanobridge junctions are incorporated into an unshunted dc 
SQUID, as shown in (b). The SQUID loop is 1.0 x 1.5 ^m, 
and is calculated to have 5 pH of inductance. 



the sample in situ and perform a second evaporation step 
to deposit 80 nm thick banks without depositing further 
material in the narrow bridge region. Atomic force mi- 
crographs of 2D and 3D nanobridges are shown in Fig. 1. 
The SQUID loop is 1.0 x 1.5 ^m, and is calculated to have 
Lioop — 5 pH of geometric inductance. Each chip con- 
tained six simultaneously fabricated nanoSQUIDs with 
bridges of varying lengths. The devices were mounted 
in a copper sample box equipped with a superconducting 
fiux-modulation coil and measured in a '^He cryostat with 
a base temperature of 265 mK. All wiring ran though 
discrete-element and distributed lossy filters to suppress 
spurious electromagnetic radiation. Four-wire transport 
measurements on all six SQUIDs were performed using 
custom electronics. 

In the inset of Fig. 2, we present a typical current- 
voltage (I-V) characteristic of a 3D nanoSQUID with 75 
nm long bridges. The I-V characteristic is hysteretic, 
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FIG. 2: (Color online) IcRn product as a function of bridge 
length for both 2D and 3D nanoSQUIDs from several chips. 
The dashed line in the main figure is the zero temperature pre- 
diction for a short metallic weak link connected to ideal phase 
reservoirs. Inset: a typical I-V curve for a 3D nanoSQUID, 
consisting of two 75 nm bridges. 



and the retrapping behavior can be attributed to ther- 
mal relaxation^. We observe that the normal state re- 
sistance Rn increases linearly vifith bridge length, with 
an offset of w 1.5 ft for 3D bridges and w 2057 for 2D 
bridges; these numbers are consistent with the normal 
state resistances of the banks. We subtract this off- 
set and plot the resulting IcRn product in Fig. 2 for 
both 2D and 3D nanoSQUIDs. For short 3D devices 
(L < 150 nm), we see that the IcRn product saturates 
at approximately 380/iF. For comparable 2D devices, 
the data are clustered around 250fj.V. For longer de- 
vices, both 2D and 3D, we see a monotonic increase 
in IcRn with bridge length, as expected for a long su- 
perconducting wire^. To interpret the data for short 
3D bridges, we can make use of the analytical result 
from Kulik and Omelyanchuk^' (KO-1), who calculate 
the current-phase relation of a weak link with L <S^ ^ 
to be I{S) = (7rA/ei?„)cos((5/2)tanh~^[sin((5/2)] where 
d is the phase difference across the junction and A is 
the superconducting gap. This formula assumes a diffu- 
sive weak link at zero temperature with a phase vari- 
ation constrained to the bridge region. The dashed 
line in Fig. 2 indicates the prediction of this relation, 
IcRn = 1.327rA/2e, which agrees well with our short 
3D bridges. For our thin films, we estimate the coherence 
length in the dirty limit to be ^ = ^JfiD/2nkBTc « 40 
nm, where D is the diffusion constant calculated from 
the measured film resistivity. Thus, our shortest devices 
are a few times the coherence length and approximate 
well a short metallic weak link connected to ideal phase 
reservoirs. 

To characterize the flux dependence of the 



nanoSQUIDs we measured I-V curves as a function 
of applied perpendicular magnetic field. From these 
data, we extracted the critical current Ic as a function 
of applied flux, and have plotted this quantity for 
nanoSQUIDS with 75, 150, and 250 nm-long bridges in 
Fig. 3. We observe that the critical current modulation 
is deeper and more nonlinear for the 3D devices. The 
modulation for the 2D devices is quite triangular in 
shape, indicative of a nearly linear current-phase re- 
lation. The data show a modulation pattern which is 
nearly symmetric about the maximum Ic- Comparing 
the measured modulation with numerical simulations 
indicates that the junction critical current asymmetry is 
typically less than 10%. Larger asymmetry would reduce 
the SQUID modulation depth and thus complicate 
systematic comparison. The modulation depth in our 
devices is also not limited by the loop inductance as 
Pl = 2LioopIc/^o ~ 0-lj where $o is the flux quantum. 
In Fig. 4, we plot the fractional modulation depth 
{Ic,max - Ic,nun)/Ic,max as a function of dcvicc length. 
For short 3D bridges (L < 150 nm), we observe a 
modulation depth approaching 70%. In longer bridges, 
the modulation depth decreases. For very short 2D 
bridges (L < 100 nm), we observe a modulation depth 
below 50%. In intermediate length 2D bridges, the 
modulation depth increases, and then decreases again 
for bridges longer than 150 nm. 

To interpret these results, let us consider a SQUID 
formed by two symmetric elements with an arbitrary cur- 
rent phase relation f{5). In the limit of small loop in- 
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FIG. 3: (Color online) Critical current as a function of applied 
magnetic flux for 75, 150, and 250 nm long nanoSQUIDs. The 
upper panel are experimental results for nanobridges with 3D 
banks and the lower panel for 2D banks. 
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FIG. 4: (Color online) Critical current modulation depth as 
a function of device length for both 3D and 2D nanobridges. 
The black dashed line is the theoretical prediction for an ideal 
weak link SQUID at zero temperature with zero loop induc- 
tance. The blue and red dashed lines are numerical calcula- 
tions for our 3D and 2D devices, respectively, including the 
effect of finite loop inductance. 



ductance, the phase difference between the two SQUID 
junctions 82 — 5i equals the applied flux 27r<I>„j/$o- The 
current / through the SQUID is the sum of the currents in 
each arm / = f{5i) + /((^s) = f{5i) + fih + 27r$^/$o). 
For a tunnel junction based SQUID, f[5) = /osin(5); 
maximizing the current equation above gives the SQUID 
critical current Ic = 2Jo| cos(7r$m/^o)| where /o is 
the critical current of each junction. The critical cur- 
rent modulates from to 2Jo and hence the modula- 
tion depth is 100%. For a linear current phase rela- 
tion, expected for a long wire, the maximum modula- 
tion depth is 25%. In the case of the KO-1 relation 
cited earlier, the maximum modulation depth is 81%. 
Thus, our short 3D bridges approach the KO-1 result, 
while the longest devices tend towards 25% modulation 



as the CPR becomes increasingly linear. Also plotted 
in Fig. 4 are the modulation depths predicted from 
our numerically computed CPRs^. These predictions 
take into account loop inductance, using the condition 
2^n = 81-82 + 27r$,„/$o + ^i/oop[/(<5i) - f{82)]l^o. 
The quantitative agreement with short 2D and 3D de- 
vices is excellent, except for the shortest 2D devices, 
which modulate less than predicted. We attribute this 
to the fact the phase gradient in short 2D bridges leaks 
into the banks of the device. Thus, for two weak links 
in close proximity it is difficult to impose a phase differ- 
ence between the two arms of the SQUID^ , and modula- 
tion is suppressed. Our calculations are for isolated junc- 
tions and a numerical calculation of the full two junction 
SQUID is required to reproduce this effect—. 

In summary, using low-inductance symmetric dc 
SQUIDs, we have demonstrated that nanobridge Joseph- 
son junctions contacted with 3D banks exhibit greater 
nonlinearity than those with 2D banks. Transport in 3D 
bridges up to several coherence lengths long is well de- 
scribed by the KO-1 formula for an ideal, short, metallic 
weak link. Our experimental results give good quanti- 
tative agreement with a numerical solution of the Us- 
adel equations for the precise geometry under study. As 
such, 3D nanobridge devices should have sufficient non- 
linearity for use in high gain amplifier and qubit circuits 
traditionally employing tunnel junctions. These devices 
could also be used for dispersive, high-speed magnetom- 
etry by embedding them in a microwave resonator—. For 
magnetometry, 3D nanoSQUIDs should have higher flux 
sensitivity than 2D devices, but perhaps with a reduced 
tolerance to in-plane magnetic fields. Further study is 
needed to determine the variation of nonlinearity and 
field tolerance with bank thickness. 
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